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Abstract 

We investigate the flux and the event rate of the supernova relic neutrino back- 
ground (SRN) at the SuperKamiokande detector for various neutrino oscillation 
models with parameters inferred from recent experimental results. A realistic model 
of neutrino emission from supernova explosions and several models of the cosmic 
star formation history are adopted in the calculation. The number flux over entire 
energy range is found to be 11 — 15 cm _2 s . We discuss the detection possibility 
of SRN at SuperKamiokande, comparing this SRN flux with other background neu- 
trinos in more detail than previous studies. Even though there is no energy window 
in which SRN is dominant, we might detect it as the distortion of the other back- 
ground event. We found in the energy range 17 — 25 MeV the expected event rate at 
SuperKamiokande 0.4 — 0.8 yr _1 . In this range, ten- year observation might enable 
us to detect SRN signal (at one sigma level) in the case of LMA solar neutrino 
solution. We also investigate event rate at SNO and KamLAND. Although we can 
find energy window, the expected event rate is rather small (0.03 yr" 1 for SNO, 0.1 
yr _1 for KamLAND). 
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1 Introduction 



A core-collapse supernova explosion produces a number of neutrinos and 99% 
of the gravitational energy is transformed to neutrinos. It is generally believed 
that the core-collapse supernova explosions have traced the star formation 
history in the universe and have emitted a great number of neutrinos, which 
should make a diffuse background. This supernova relic neutrino (SRN) back- 
ground is one of the targets of the currently working large neutrino detectors, 
Super Kamiokande (SK) and Sudbury Neutrino Observatory (SNO). Compar- 
ing the predicted SRN spectrum with the observations by these detectors 
provides us potentially valuable information on the nature of neutrinos as 
well as the star formation history in the universe. This SRN background has 
been discussed in a number of previous papers [1]-[11]. The work after Totani 
et al. [8] takes into account the realistic star formation history inferred from 
various observations and theoretical modeling of galaxy formation, to calcu- 
late the SRN flux. Totani et al.[8] calculated the energy flux of SRN at SK 
detector and compared it with neutrinos emitted by other sources (solar, at- 
mospheric, and reactor neutrinos and so on). Then they concluded that the 
visible event rate of SRN at SK is 1.2 yr" 1 in the energy range from 15 to 40 
MeV. On the other hand, Kaplinghat et al. [11] calculated the upper limit of 
SRN and discussed the possibility of the detection at SK detector comparing 
to the other backgrounds. Their result is that at the energy range from 15 
to 40 MeV, where no significant background had been considered to exist in 
previous studies, there is a huge background of the invisible muon decay and 
the detection of SRN is difficult. They also discussed the effects of neutrino 
oscillation and briefly mentioned that if the SRN flux is in the vicinity of their 
upper bound and all three flavors are maximally mixed, it may be detectable 
as a distortion of the expected muon background. 

In this paper we calculate the SRN flux and the event rate at SK, and discuss 
about the detectability of SRN, with the following new aspects compared with 
previous studies: 1) realistic neutrino oscillation parameters are incorporated 
based on the recent solar and atmospheric neutrino experiments, 2) a realistic 
neutrino spectrum from one supernova explosion is used, which is obtained 
from a numerical simulation by the Lawrence Livermore group, and 3) we 
have examined other contaminating background events against the detection 
of SRN, in more detail than previous studies. 

Recent experiments of SK and SNO support neutrino oscillation by solar [12]- 
[14] and atmospheric neutrino [15] data. Although some previous studies (e.g. 
[11]) mentioned the possibility of neutrino oscillation, no quantitative calcula- 
tion of SRN has been made incorporating the realistic oscillation parameters. 
Here we consider four neutrino oscillation models, which satisfy the solar and 
atmospheric neutrino data, and investigate the dependence of the positron 
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spectra at SK on the oscillation models. If neutrino oscillation occurs, z/ MiT 's 
are converted into z/ e 's which are mainly detected at SK detector. Because 
z/ MjT 's interact with matter only through the neutral-current reactions in su- 
pernovae, they are weakly coupled with matter compared to V e . Thus the 
neutrino sphere of z/ MiT 's is deeper in the core than z/ e 's and their temperatures 
are higher than z/ e 's. Therefore neutrino oscillation enhances the mean P e 's 
energy and enhances event rate at SK detector. 

In all of the past studies, neutrinos emitted from a supernova are assumed 
to obey the Fermi-Dirac distribution with zero chemical potential. However, 
since neutrinos are not in the thermal equilibrium states in supernovae, the 
real spectrum should be different from the pure Fermi-Dirac distribution. Thus 
we use in this paper a realistic supernova model established by the Lawrence 
Livermore group [16] which in fact shows clear difference from the Fermi-Dirac 
distribution [17]. 

There are several background events which hinder the detection of SRN. These 
includes atmospheric and solar neutrinos, anti-neutrinos from nuclear reactors, 
and decay electrons from invisible muons. We should find the energy region 
which is not contaminated by these background events and then calculate the 
detectable event rate of SRN. By careful examination of these events, we found 
that there is a narrow energy window at which the detection of the SRN event 
might be possible. 

In addition we use three models of the cosmic star formation history, i.e., the 
evolution of star formation rate (SFR) density of high redshift galaxies (e.g., 
[18]), which is inferred from several star formation indicators such as rest-frame 
UV luminosity, Ha lines, or dust emission in far-infrared or submillimeter 
wavebands. It should be noted that SFR models contain various uncertainties, 
especially at high redshift regions, {e.g., the faint ends of luminosity functions 
have not been well established at high redshifts and the uncertainties in dust 
extinction are large at all redshifts [19].) This is the reason we applied several 
SFR evolution models described below. 

Throughout this paper, we consider only electron antineutrinos (P e ' s ) from 
collapse-driven supernovae, because P e 's are most easily detected in a water 
Cherenkov detector like SK. 

This paper is organized as follows: In Section 2, we illustrate the neutrino oscil- 
lation models, the supernova model, and the supernova rate models considered 
in this paper and discuss about these models. Formulations for calculation of 
flux and event rate at SK are given in Section 3. In Section 4, we show the 
calculated flux and event rate. Detailed discussions especially on other back- 
ground events are presented in Section 5. 
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2 Models of Neutrinos and Supernova Rate 

2.1 Supernova Model 

We use a realistic model of a collapse-driven supernova calculated by the 
Lawrence Livermore group [16]. We show in Fig. 1 the time-integrated spec- 
trum for UeS (e.g., see Ref. [17] for detail). We also show the Fermi-Dirac (FD) 
distribution with zero chemical potential for comparison. Comparing with the 
FD model used in the other studies, the deficit of both low- and high-energy 
neutrinos can be seen. 

In this paper, we assume that this numerical model represents all of the past 
supernovae, although in this model the progenitor of the supernova has been 
assumed to have the mass of ~ 2OM . It is clearly an oversimplification, but we 
note that the mean mass of progenitor stars of type II supernovae (above 8M ) 
weighted by number is about ~ 15M Q when a typical initial mass function is 
applied [3]. 

Another assumption made in the simulation is the isotropic radiation of neu- 
trinos. All supernova progenitors are rotating, and it may have significant 
effect on the degree of isotropy of neutrino emission [20], although the effect 
is very difficult to estimate quantitatively. However, since our interest is in 
the SRN background, which is the sum of all past supernova neutrinos, the 
rotation effect is expected to be small. 

2.2 Neutrino Oscillation Models 

We adopt the four parameter sets for the neutrino mixing (see Table 1) and 
the normal mass hierarchy. These parameter sets are introduced to explain the 
observations of the solar and the atmospheric neutrinos [12] -[15]. The neutrino 
spectra emitted by each supernova are calculated numerically by Takahashi 
et al. [21] for the same parameter sets assuming the normal mass hierarchy. 
(Supernova model they have used is the same one discussed in the previous 
subsection.) We use their calculated u e spectrum as that from any supernova. 

In Table 1, "LMA" and "SMA" indicate MSW solution of the solar neutrino 
problem (see Ref. [22], for the review of MSW effect). Recent SK and SNO 
observations [12,14] show that LMA solution are more favorable but SMA so- 
lution cannot be rejected at more than 3<r level (see Ref. [23] and references 
therein). Therefore we try the SMA solution as well. The suffixes "-L" and 
"-S" indicate whether #13 is large or small. A large (small) #13 means "higher 
resonance" is adiabatic (nonadiabatic). (There are two resonance points where 
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neutrino mass eigenstates may flip in the supernova matter, and we call them 
lower and higher resonance points. The lower resonance point is at lower den- 
sity, and the higher resonance point at higher density.) The adiabatic higher 
resonance enhances the energy of the electron neutrinos and enhances the 
event rate of v e scattering, but does not affect the electron anti-neutrinos [24]. 
In this paper, since we deal with only the reaction of z/ e 's, this "-L" and "-S" 
hardly influence the result. 

Now, we semi-qualitatively present a simple illustration of the expected spec- 
tral shape of P e 's for LMA and SMA models. We can naively deal with anti- 
neutrino oscillation effect as vacuum oscillation, since z/ e 's are not affected 
by the resonance. Further we take two generation formalism, for simplicity 
(this approximation is justified when # 13 is sufficiently small). The conversion 
probability from to u e (and its inverse) averaged over distance is 

Pfa -> Ve) = P(p e -> P M ) = - sin 2 2012. (1) 



This probability depends only on 9 12 . Then, P{v^ — > 9 e ) = P{y e — > v^) = 
0.44(LMA), 2.5 x 10~ 3 (SMA). Using these conversion probabilities we can 
present V e flux as 



[1 - P{u e -> V,)\Fl + P(P, -> u e )Fl 

0.57F P ° e + 0.44F p ° m (LMA) 
K (SMA) ' 



(2) 
(3) 



where F® means neutrino flux for no oscillation model. This equation shows 
that in LMA models the v e flux at higher energy region is enhanced, on the 
other hand that at lower energy region is suppressed, because the original 
mean P^s energy is higher than that of z/ e 's as mentioned in Section 1. 

However, in the case of the inverted hierarchy, this situation changes dra- 
matically although we do not give calculation for this case. In this case, the 
higher resonance point is in the anti-neutrino sector [24], and hence resonant 
conversion of z/ e 's and is possible. The conversion probability would be 

dependent on the neutrino energy. Another possibility of such resonant con- 
version effect is the spin-flip oscillation of z/ e 's by a flavor- changing neutrino 
magnetic moment [25]. 
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2.3 Supernova Rate 



A number of studies have modeled the expected evolution of the cosmic SFR 
with redshift. The cosmic SFR can now be traced to z ~ 4 observationally, 
although some details remain controversial. We use three SFRs per unit co- 
moving volume which are also used in Ref. [26], where all of them are derived 
following the same method employed by Madau et al. [18] first. They use the 
intergalactic absorption to identify high redshift galaxies in broadband mul- 
ticolor survey, and convert the UV luminosity to SFR and the metal ejection 
rate assuming a typical initial mass function (IMF) of stars. The supernova 
rate is expected to be proportional to SFR since the lifetime of progenitors of 
core-collapse supernovae is much shorter than the cosmological time scale. We 
show in Fig. 2 the supernova rates corresponding to these three SFR models 
that are explained below. 

In the Einstein-de Sitter (EdS) universe (Q m = 1.0, = 0.0), the first (here- 
after SF1) is taken from Ref. [27]: 

W«) = 0.3 ^g^ a yr--Mpc-', (4) 



where = Ho/65 km s -1 Mpc -1 . This SFR increases rapidly between z = 
and z — 1, peaks between z — 1 and z — 2, and gently declines at higher 
redshifts. This includes an upward correction for dust reddening of A 1500 = 1.2 
mag. The original SFR (before correction) should be interpreted as lower limits 
to the real values, since a significant fraction of UV light from young stars 
could be absorbed and re-emitted in far-infrared band. In the second model, 
the SFR remains instead roughly constant at z > 2 (SF2) [19], as: 

R SF2 (z) = 0.15/i 65 ^ 3 ; 4 j 00 M Q yr^Mpc- 3 . (5) 
exp(3.4z) + 22 



We should consider SF2 in addition to SF1, because of the uncertainties as- 
sociated with the incompleteness of the data sets and the amount of dust 
extinction at early epochs. The third SFR (SF3) [26], 



„ . , , exp(3.05z - 0.4) w , , 



(6) 



represents even more star formation at early epochs. This SFR is based on 
the studies which suggest that the evolution of the SFR up to z ~ 1 may have 
been overestimated [28], while the rates at high-z may have been severely 
underestimated due to large amounts of dust extinction [29]. 
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These SFR evolutions should change when a different cosmological model is 
assumed, since they are based on the observed data of high- z galaxies. The 
correction to other cosmological models can be written as (e.g., [30]): 



RSF{Z] &m, ^65) — h$5 



(l + n m z)(l + z) 2 -tt x (2z + z 2 ) 



(1 + ,2)3/2 

x R SF (z; 1,0,1). 



(7) 



To obtain the supernova rate (Rsn), we multiply the SFRs by the coefficient 

r/^r' -0.0122V. < 8 » 

J0.01 dmm(f)(m) 



where <f>(m) is the Salpeter IMF (0(m) oc m~ 2 - 35 ) and m is the stellar mass 
in solar units [26]. (Here, we assume that all stars whose mass is greater than 
8M explode as core-collapse supernovae.) This resulting rates agree with 
the local observed value. We label these models "SN1", "SN2", and "SN3", 
respectively. 



3 Formulation of Flux Calculations 



In section 2, we calculated the number of emitted z/ e 's from a supernova per 
unit energy q, dN u /dq, and supernova rate, Rsn- We remark that Rsn(z) is 
supernova rate per comoving volume, and hence we should multiply the factor 
(1 + z) 3 to obtain the rate per physical volume at that time. The present 
number density of P e 's, whose energy is in the interval of the q ~ q + dq, 
emitted in the interval of the redshift z ~ z + dz is given by 



dn v {q) = Rsn{z){1 + zf^dz dNv ^\ + z)q \ l + z)dq{\ + z)~ 3 (9) 

dz dq 

„ / > dt dNJ(l + z)q) , 
= R SN (z)—dz )q, (l + z)dq, (10) 

where the factor (1 + z)~ 3 comes from the expansion of the universe. The 
Friedmann equation gives the relation between t and z as follows: 

dz 



— = -h q (i + ^(1 + n mZ )(i + zf - n x (2z + z*). (11) 
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We now obtain the differential number flux of SRN, dF l/ (q)/dq, using the 
relation dF u (q)/dq = c[dn u (q) / dq}: 



dF v c 



dq H Q 




(12) 



where we assume that gravitational collapses began at the redshift parameter 



It can be seen that, from equations (7) and (12), the flux does not depend 
on the cosmological parameters such as fl m , Q\, and H , while in Ref. [8], the 
flux depend on these parameters. The illustration of this difference is as fol- 
lows: Totani et al. used supernova rate evolution derived from their theoretical 
model of galaxy evolution which reproduces various properties of present-day 
galaxies. On the other hand, our supernova rate is based on the observational 
estimate of luminosity densities of high- z galaxies, and the dependence of the 
cosmological volume element on the parameters such as fl m and Q\ is cancelled 
out. Therefore, our supernova rate depends on the cosmological parameters 
while the neutrino flux does not depend on them. 



4 Results 

4-1 Calculation of SRN Flux 

We calculate the flux of the SRN using the formula eq. (12) for various models. 
In Fig. 3, the flux for the three supernova rate models are shown assuming no 
oscillation model. The fluxes of these models are almost the same above ~ 8 
MeV and the models with more SFR at early epochs have a higher peak at 
lower neutrino energy. These properties come from the following effects: The 
energy of neutrinos which were emitted at redshift z is reduced by a factor 
(l+z)" 1 when we observe. Then at observation, the high energy tail (> 8 MeV) 
is mainly contributed by low redshift supernovae. Since at low redshift three 
supernova rate models are almost the same, the fluxes at E Pe > 8 MeV are not 
much different. Similarly, because high redshift supernova neutrinos contribute 
more at low energy region when we observe, we can see model dependence at 
low energy clearly. 

In Fig. 4, we show the flux of various neutrino oscillation models assuming 
SN1 model. In SMA and no oscillation models we can see higher peak around 



in ax 



5. 
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~ 5 MeV than in LMA models, because in LMA models the low energy z/ e 's 
are deficit due to the conversion into z^ )T 's. But above ~ 10 MeV, in LMA 
models we can see more flux, because z/^ T 's which have higher mean energy 
at production have more changed into z/ e 's than in the case of SMA or no 
oscillation. (See also qualitative explanation in Section 2.2.) 

Integrated flux over the entire neutrino energy range is shown in Table 2. 
4-2 Event Rate at SuperKamiokande Detector 

The SK detector is a water Cerenkov detector whose fiducial mass is 22,500 
ton. The detector efficiency is 100% for electrons (positrons) whose energy is 
above 5 MeV, and 50% at 4.2 MeV. We only consider the reaction v e + p — > 
e + + n, because the cross section [9.52 x 10~ 44 (i? e /MeV)(p e /MeV)cm 2 [31]] is 
much larger than the other reactions. 

In Fig. 5, the event rate at SK of the three supernova rate models are shown 
assuming no oscillation model. Because of the detector threshold, we see only 
the positrons whose energies are above ~ 5 MeV, so that the differences be- 
tween models are small as shown also in Fig. 3 where the model dependence is 
also small for neutrinos above ~ 7 MeV. Then, although our three supernova 
rate models have very different properties at high redshift regions, our results 
are not influenced by the behavior of SFRs at those high z regions. (These 
results are hardly changed when we assume the other LMA or SMA models, 
instead of the no oscillation model.) 

In Fig. 6, we show the event rate of various neutrino oscillation models assum- 
ing SN1 model. We can see clear difference of the LMA model from SMA or 
no oscillation models, especially at the high energy tail. This property results 
from the flux dependence on oscillation models (see also Fig. 4). 

As a result, if we can detect SRN events above ~ 10 MeV, we can discriminate 
the LMA model from the SMA or no oscillation models in any supernova rate 
models and any sets of cosmological parameters. 



5 Discussion 

5. 1 Background Events against the Detection 

We discuss in this subsection about neutrinos from other sources which may 
become an obstacle to the SRN detection. They are atmospheric and solar 
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neutrinos, anti-neutrinos from nuclear reactors, and decay electrons from in- 
visible muons. We show in Fig. 7 the number flux of SRN and these background 
events, and we discuss the each background event below. 

The flux of the atmospheric neutrinos is usually calculated using Monte Carlo 
method including various relevant effects (flux of primary cosmic rays, solar 
modulation, geomagnetic field, interaction of cosmic rays in the air, and so on), 
and in that simulation one-dimensional approximation is used, i.e., after the 
interaction of primary cosmic ray particles with air nuclei, all the particles are 
assumed to be moving along the line of the momentum vector of the primary 
cosmic ray particles [32]. ( Recently preliminary results of three-dimensional 
flux calculations have been reported [33].) There are many authors who cal- 
culated the atmospheric neutrino flux (see, e.g., [34] for a recent result). We 
use in this paper the flux calculated by Gaisser et al. [35] and Barr et al. [36]. 
Although their calculations are rather old, the fluxes of low energy neutrinos, 
in which we are interested, are also calculated, while most of other papers 
show only higher energy region (> 1 GeV). 

Solar neutrino flux is dominant at energy range 10-20 MeV. We use the flux 
predicted by the standard solar model (SSM) in Fig. 7 [37]. Since the solar 
neutrinos are not z/ e 's but z/ e 's, the cross section for solar neutrinos is about 
two order smaller than that for z/ e 's. Furthermore recoil electrons scattered by 
solar neutrinos strongly concentrate to the opposite direction of the Sun, in 
contrast to the isotropic distribution of v e events. Therefore the solar neutrino 
is an avoidable background. In fact, we can demonstrate how we can avoid 
the solar neutrino background, as follows. The solar neutrino event number is 
about 5 x 10 3 yr~ 1 (22.5kton) _1 [12], about 3 or 4 orders of magnitude larger 
than our expected event rate ~ 1 yr _1 (22.5kton) _1 . The recoil electrons scatter 
obeying the Gaussian of one-sigma error o ~ 25° for 10 MeV electrons [38]. 
Then, 99.98% of solar neutrino events within ~ 3.7cr (corresponding to ~ 90°) 
can be avoided. Therefore, restricting our discussion on the nearby side of the 
hemisphere from the Sun, we can ignore the solar neutrino events. 

The third background which we must consider is anti-neutrinos from nuclear 
reactors. In each nuclear reactor, almost all the power comes from the fissions 
of the four isotopes, 235 U (~ 75%), 238 U (~ 7%), 239 Pu (~ 15%), and 241 Pu 
(~ 3%) [39]. Each isotope produces a unique electron anti- neutrino spectrum 
through the decay of its fission fragments and their daughters. The v e spectrum 
from 235 U, 239 Pu, and 241 Pu can be derived using the semi-empirical formula 
with which we fit data of detected /5-spectrum from fission by thermal neutrons 
[40]. ( 238 U undergoes only fast neutron fission and hence electron spectrum 
from 238 U cannot by measured by this kind of experiment.) Above 7 MeV, 
the number of (3 counts drops dramatically and fitting error becomes large. 
In addition, with this method, as we determine the maximum /3 energy and 
derives the energy distribution below that energy, it is difficult to estimate 
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the errors on the high energy range [41]. While the v e spectra in Ref. [40] 
are given as tables, we use for simplicity somewhat less accurate analytical 
approximation given in Ref. [42]. As a normalization factor we use energy- 
integrated V e flux at Kamioka, 1.34 x 10 6 cm~ 2 s _1 , which are the summation 
of the flux from various nuclear reactors in Japan and Korea [39]. This fit is 
not valid above ~ 8 MeV. However, there is an estimation that we would get 
few (10 or less) events per year above 10 MeV from reactors[41]. 

With three backgrounds we discussed above, we expect the energy window of 
SRN events opening from 10 MeV to 27 MeV. However, according to Kap- 
lighat et al. [11] electrons or positrons from invisible muons are the largest 
background in the energy window from 19 to 35 MeV. This invisible muon 
event is illustrated as follows. The atmospheric neutrinos produce muons by 
interaction with the nucleons (both free and bound) in the fiducial volume. If 
these muons are produced with energies below Cherenkov radiation threshold 
(kinetic energy less than 53 MeV), then they will not be detected ("invisible 
muons"), but their decay-produced electrons and positrons will be. Since the 
muon decay signal will mimic the u e p — > ne + process in SK, it is difficult 
to distinguish SRN from these events. The energy spectrum of this invisible 
muon events is obtained by the stopped muon decay spectrum 



(Michel spectrum [43]), where G F is the Fermi constant, the muon mass, 
and E e the electron (positron) energy. This equation is valid for E e < 
From the observation of Kamiokande II detector the estimated event rate 
from these muon decays is around unity for 0.58 kton yr exposure and this 
is forming the principal source of background after the various cuts had been 
implemented [44]. Kaplinghat et al. [11] used this value corrected for SK by 
multiplying the volume ratio and concluded that it is impossible to detect SRN 
unless their upper limit realizes and maximum neutrino oscillation occurs. This 
event rate is now inferred directly from the new data of the SK, which is about 
100 per 1,258 days per 22.5 kton fiducial volume above 18 MeV [41,45]. This 
value corresponds to about 20% of the Kamiokande II data. In Fig. 8 we show 
SRN event rate compared to the invisible muon events. From this figure we 
conclude that SRN events can be seen only below about 12 MeV. 

In practice, there is another serious background, i.e., spallation products in- 
duced by cosmic ray muons. Ultra high energy cosmic ray muons spall oxygens 
in the detector, and radioactive decay processes of these spalled nuclei occur. 
The event rate of the spallation background is several hundred per day per 22.5 
kton. Although most of them can be rejected by the information of preceding 
muons, even a small fraction can not be. (Roughly, this spallation products 
produces about 200,000 events per a year. Because expected SRN event rate is 




(13) 
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less than 1 per a year, we should reject all these 200,000 events. For future de- 
tectors this problem is also quite difficult to solve.) Then, this makes a serious 
background at the energy range below the maximum energy of beta spectrum 
of spallation products, 16 MeV [45,46]. From these discussions presented in 
this subsection we conclude that there is no energy window of SRN. 

5.2 Calculating the detectable event rate at SK 

In the previous subsection we have found that there is no energy region where 
SRN is dominant. However, we can detect the SRN events by subtracting the 
other background events from total detected events. We consider the energy 
range 17 < (T e /MeV) < 25, where T e is positron kinetic energy. This range 
corresponds to 19 < (E Pe /MeV) < 27 by the simple relation, E Pe = T e + 
1.8MeV. We find two advantages in using this energy region. First, SRN event 
rate is rather large, and second, the background (invisible muon) event rate 
is fairly well known by SK observation. We show in Table 3 the SRN event 
rate at SK in this energy range, i.e., 0.4 — 0.8 yr _1 . In contrast, the event 
rate of the invisible muon over the same energy range is 3.4 yr _1 . When 
SRN event rate is larger than the statistical error of background event rate, 
we can conclude that the SRN is detectable as a distortion of the expected 
invisible muon background event. Unfortunately, only one year observation 
does not provide any useful information about SRN. However, we can expect 
that ten- year observation provides several statistically meaningful results. The 
statistical error of invisible muon events in ten years is \^3A = 5.8, which is 
smaller than the event rate of LMA models and is larger than that of SMA 
and no oscillation models. Then we conclude that these neutrino oscillation 
models can be distinguished by the observation of the event rate of invisible 
muon events. (If there is a discrepancy from expected event rate, this is due 
to SRN events and LMA models are favored.) 

In future, it is expected that next generation of water Cherenkov detectors 
have much larger volume than that of SK. For example HyperKamiokande 
project is now under consideration. HyperKamiokande detector is planed to 
be a water Cherenkov detector whose mass is about 1,000,000 tons (about 20 
times larger than SK), and its location is near SK detector. We expect that the 
SRN event becomes about 10 per one year for this detector, and statistically 
sufficient discussion of SRN is possible even using only one year data. 

5.3 SNO and KamLAND detectors 

In this subsection we discuss SRN detectability with SNO [47] and KamLAND 
[48] . An advantage of these detectors is that we are able to identify v e events 
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using delayed coincidence signals. For this reason, we can remove other back- 
grounds from non-z/ e origin (solar neutrinos, invisible muon decay products, 
and spallation products). 

At SNO detector, used reaction for v e detection is 

V e + d — > e + + n + n. (14) 



Then these neutrons react with surrounding nuclei through 



n + d — > 3 H + 7 (E y = 6.3MeV, efficiency : 24%), (15) 
n + 35 C1 — ► 36 C1 + 7 (£ 7 = 8.6MeV, efficiency : 83%), (16) 

where NaCl is added to efficiently capture neutrons. This delayed signal with 
the preceding Cherenkov radiation from e + in eq. (14) shows that the detected 
neutrino is v e . Using this criterion we can reject backgrounds from non-z/ e 
origin. These are solar neutrinos, invisible muon decay products, and spallation 
products, which are great obstacle for SRN detection at SK. Thus, we can 
find the energy window below T e ~ 23MeV. Actual calculation shows that the 
expected event rate at SNO is ~ 0.03 yr -1 . (We used cross section calculated 
in Ref. [49].) 

We can use a similar criterion at KamLAND. P e 's are detected through below 
reactions: 

v e +P — >e + +n, (17) 
n + p — > d + 7 (E y = 2.2MeV). (18) 



The energy window for SRN is from 10 MeV to 25 MeV. (Below 10 MeV 
reactor P e 's are large.) In that range the calculated event rate is ~ 0.1 yr -1 . 

Unfortunately these values (0.03 yr -1 for SNO; 0.1 yr -1 for KamLAND) are 
quite small, since the fiducial volume of these detectors (1 kton) is much 
smaller than that of SK. However, the future same kind of detector of larger 
volume, if ever built, might detect SRN. 
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Neutrino Energy [MeV] 

Fig. 1. Energy spectrum of v e of the numerical supernova model used in this paper. 
The dashed line is the Fermi-Dirac fits which have the same luminosity with the 
numerical model. The chemical potential is set to zero for the FD distribution. 
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Fig. 2. Supernova rate evolution on the cosmological time scale. These lines are for 
a A-dominated cosmology (f2 m = 0.3, = 0.7). The Hubble constant is taken to 
be 70 km s _1 Mpc -1 . 
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Fig. 3. Number flux of P e 's for the three supernova rate models, assuming "no 
oscillation" case. 
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Fig. 4. Number flux of P e 's for the neutrino oscillation models. In this figure "SN1" 
model is used for supernova rate evolution. 
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Fig. 5. Event rate of z/ e 's at SK for the three supernova rate models, assuming "no 
oscillation" case. 
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Fig. 6. Event rate of P e 's at SK for the neutrino oscillation models. In this figure 
"SN1" model is used for supernova rate evolution. 
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Fig. 7. Number flux of SRN compared to other background neutrinos. No oscillation 
and SN1 model are assumed for SRN flux. 
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Fig. 8. Event rate at SK detector of SRN and invisible /j, decay products. Two 
oscillation models are shown (no oscillation and LMA-L) assuming SN1 model for 
supernova rate evolution. 
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Table 1 



Sets of mixing paremeter for calculation 



model 


sin 2 26»i2 


sin 2 26*23 


sin 2 26*13 


Am 2 2 (eV 2 ) 


Am 2 3 (eV 2 ) 


Vq problem 


LMA-L 


0.87 


1.0 


0.043 


7.0 x 10~ 5 


3.2 x 10- 3 


LMA 


LMA-S 


0.87 


1.0 


1.0 x 10~ 6 


7.0 x 10~ 5 


3.2 x 10~ 3 


LMA 


SMA-L 


5.0 x 10~ 3 


1.0 


0.043 


6.0 x 10~ 6 


3.2 x 10^3 


SMA 


SMA-S 


5.0 x 10~ 3 


1.0 


1.0 x 10~ 6 


6.0 x 10~ 6 


3.2 x 10^3 


SMA 



Table 2 

Number flux of SRN. Each entry is the number flux in unit cm^ 1 s _1 for each 
supernova rate and oscillation models. These fluxes are integrated over whole energy 
range. 



models 


LMA-L 


LMA-S 


SMA-L 


SMA-S 


no oscillation 


SN1 


11.2 


11.3 


12.3 


12.2 


11.9 


SN2 


12.2 


12.2 


13.4 


13.2 


12.9 


SN3 


13.8 


13.8 


15.1 


14.9 


14.6 



Table 3 

SRN event rate at SK detector. Each entry is the event rate in unit yr~ x for each 
supernova rate and oscillation models. Integrated energy range is from 17 to 25 
MeV. 



models 


LMA-L 


LMA-S 


SMA-L 


SMA-S 


no oscillation 


SN1 


0.73 


0.72 


0.46 


0.45 


0.44 


SN2 


0.69 


0.68 


0.44 


0.43 


0.42 


SN3 


0.76 


0.75 


0.50 


0.49 


0.48 



25 



